Phonon anomalies due to strong electronic correlations in layered organic metals 
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We show how the coupling between the phonons and electrons in a strongly correlated metal 
can result in phonon frequencies which have a non-monotonic temperature dependence. Dynamical 
mean-field theory is used to study the Hubbard-Holstein model that describes the k-(BEDT-TTF)2X 
family of superconducting molecular crystals. The crossover with increasing temperature from a 
Fermi liquid to a bad metal produces phonon anomalies that are consistent with recent Raman 
scattering and acoustic experiments. 
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An important problem concerning strongly correlated 
metals such as heavy fermions, cuprates and organic 
superconductors is understanding the interplay of the 
strong interactions between electrons and the interac- 
tions between the electrons and phonons. This is partic- 
ularly relevant to understanding the question of whether 
the phonons play any role in superconductivity ffl. In 
this Letter we show how strong electronic correlations 
can lead to phonon frequencies varying with temperature. 
Although, our calculations focus on explaining recent ex- 
periments on a particular family of organic superconduc- 
tors, the physics involved is relevant to other strongly 
correlated metals in which there is a significant redistri- 
bution of the electronic spectral weight as the tempera- 
ture is varied. 

The quasi-two-dimensional organic superconductors, 
k-(BEDT-TTF)2X, are strongly correlated electron sys- 
tems Recent Raman scattering experiments 
find that in the metallic state the frequency of certain 
phonons associated with the BEDT-TTF molecules have 
a non-monotonic temperature dependence below 200 K. 
Acoustic experiments || find that around 40 K there is 
a significant softening of the speed of longitudinal sound 
propogating perpendicular to the layers. This softening, 
of the order of a few per cent, is more than an order of 
magnitude larger than the softening associated with the 
superconducting transition. We will show that such a 
temperature dependence can arise due to the destruction 
of Fermi liquid quasiparticles that occurs above the co- 
herence temperature, T* P,^0|,p[. Anomalies in acoustic 
phonons in the heavy fermion UPt3 have also been seen 
at temperatures of the order of T* [|llj . 

The simplest possible strongly correlated electron 
model for the k-(BEDT-TTF) 2 X family is a single band 
Hubbard model on an anisotropic triangular lattice at 
half-filling ||0]. We have found || that many of the 
transport properties of the metallic phase are consis- 
tent with the predictions of dynamical mean-field theory 
(DMFT) H which captures exactly the dynamical fluc- 
tuations at each lattice site, but neglects all non-local 
spatial correlations. In order to understand the effect of 



electronic correlations on the molecular phonon modes 
we have considered the relevant Holstein-Hubbard model 
where the phonon amplitude couples linearly to the local 
charge density. The phonon self-energy is proportional to 
the electronic density-density correlation function which 
we calculate using DMFT. When the electron-electron 
interactions are sufficiently large we find that phonon 
frequencies can have a non-monotonic temperature de- 
pendence. This is related to the crossover from a Fermi 
liquid to a bad metal that occurs with Increasing tem- 
perature and has a significant effect on the electronic 
transport properties Q . The actual form of the temper- 
ature dependence of the shift in phonon frequency varies 
significantly with the phonon frequency. 

We consider the Hubbard-Holstein Hamiltonian 

H = ti J^(4 a c j(T + h.c.) + t 2 J2(4a c k* + h.c.) 

i ia 
ia i 

where the electronic part describes electrons on the anti- 
bonding orbitals of each dimer of BEDT-TTF molecules 
which are located on an anisotropic triangular lattice. t\ 
and ti are the nearest and next-nearest-neighbour hop- 
pings. U is the Coulomb repulsion between two electrons 
on the same site and y, is the chemical potential. The 
operator, cj, creates an electron on the anti-bonding or- 
bital of the dimer. The operator, aj, creates a phonon at 
site i, which describes a molecular vibration of frequency 
luq ■ We consider only the in-phase vibrations of the two 
phonon modes associated with the dimer; these can be 
activated by Raman scattering (see below), g is the cou- 
pling between the Raman active phonon and the electron 
density on a dimer. 

We focus on the parameter regime where the electron- 
electron interaction is dominant and we are well away 
from any instability (superconducting or charge-density- 
wave) due to the electron-phonon coupling. We can then 
decouple the set of Dyson's equations in the electron- 
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phonon problem, so that the electron self-energy contains 
the electron-electron scattering mechanism only; any ef- 
fects coming from the interaction of the electrons with 
phonons on the electron propagator are neglected. The 
electron Green's function is given by |9] 



G a {k,iw n ) = 



1 



6k - E(Kt>n) 



(2) 



where tu n = (2n+l)irT is a Matsubara fermion frequency 
for temperature T. — t\ cos(k x ) + t% cos{k x + k y ) is 
the dispersion relation for the anisotropic triangular lat- 
tice. Yi{iuj n ) is the momentum independent self-energy 
computed within DMFT from the associated Anderson 
impurity problem, using the iterative perturbation the- 
ory. Details can be found elsewhere J||| . 

The phonon problem is solved separately through the 
associated Dyson equation 



.g 2 ^n(q,w)/2 



(3) 



where n(q, lu) is the electronic density-density correlation 
function (or polarization) which includes the full effect of 
the electron-electron interactions at the level of DMFT. 
For n(q, ui) we take the particle-hole bubble, which, in 
terms of Matsubara frequencies is given by 

n(q,iv n )=T ^2 G a (k,iuj n )G !7 (k + ci,iuj n + ii' n ) (4) 



and G a (k,iuj n ) is the Green's function of the electrons 
obtained from the solution of (||). The temperature 
dependence of the polarization predominantly comes 
from the temperature dependence of the individual one- 
electron Greens functions. Fig. |l| shows how the Fermi 
liquid quasiparticle peak in the spectral density of states 
is strongly temperature dependent. We show results for 
the case, t\ = t% = t j^J, which corresponds to a trian- 
gular lattice. 

Due to the electron-phonon interactions, the phonon 
frequency is shifted from its bare value loq and has a fi- 
nite lifetime. This shift can be obtained from the poles 



of 



which we denote by (D q 



iTq. For weak 



electron-phonon coupling the frequency shift is 



= 2-n*(q,a, ) 
uj 4 



(5) 



where n fl (q, u>q) denotes the real part of the polariza- 
tion. The phonon damping is proportional to the imag- 
inary part of the polarization. Fig. ^| shows the tem- 
perature dependence of the real part of the polarization, 
II- R (q = 0, wo); for different bare phonon frequencies, u>o- 
Our DMFT calculations show that for frequencies close 
to U /2 and sufficiently strong interactions the phonon 
frequency can have a non-monotonic temperature depen- 
dence near the coherence temperature, T* . From Fig. |^, 



we find that T* w O.lt ~ 100-fT, for t ~ 0.1 eV. However, 
this behaviour disappears gradually as U decreases and 
becomes smaller than the bandwidth, W = 4.5i. This 
is clearly seen for U = 2t. The effective mass enhance- 
ment for U = 5.5t is m*/m w 3.8; this is in the range 
of the effective masses found in k-(BEDT-TTF) 2 X [|. 
Lin et al. [Q have measured an anomalous softening be- 
low about 100 K of the Raman frequency shifts in the 
phonons, z^g = 505 cm" 1 and V6 (B 3g ) — 890 cm" 1 for 
k-(BEDT-TTF) 2 Cu(SCN) 2 (see Fig. 5 in ref. §). The 
higher frequency mode, v% = 1478 cm" 1 , also exhibits a 
strong (but monotonic) temperature dependence. Table 
[| gives the magnitude of the total temperature depen- 
dence between about 10 and 300 K. Other phonons do 
not exhibit such a strong temperature dependence [fl3|| . 
Fig. H and equation (^) implies that for a coupling g ~ t 
the temperature dependence can be as large as 5 % for 
ujq = [7/2 and of the order of 1 % for larger or smaller 
frequencies. The values of g deduced below (see also Ta- 
ble 1) are consistent with this if t ~ 0.1 eV which is also 
reasonable |J. However, if g ~ t ~ uj this raises ques- 
tions about vertex corrections to ^ and the contribution 
of electron phonon coupling to the electronic self energy. 
Yet the effects shown in Fig. || are at most a few per cent 
and so we consider that effects that are higher order in g 
will be very small. 

We find even stronger effects for low-frequency 
phonons. In Fig. 3 the real part of the phonon self- 
energy |Q is plotted for q = and lu = O.Obt. For 
strong correlations (U > W), a dip in the real part of 
the self-energy appears at the coherence temperature, 
T*. For decreasing values of U, the position of the dip 
moves to higher temperatures (because T* increases) and 
the dip becomes smaller. These results could be rele- 
vant to understanding recent acoustic experiments [0. 
The velocity of ultrasonic waves which are propagating 
perpendicular to the layers was found to have a non- 
monotonic temperature dependence. These waves have 
frequencies of 100 MHz and velocities of about 2000 cm/s. 
The velocity versus temperature shows a broad dip of 
a few per cent around 40 K. This softening becomes 
larger as the pressure is decreased and is about three 
times larger for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br than k- 
(BEDT-TTF) 2 Cu(SCN) 2 . Decreasing the pressure or 
changing the anion from Cu(SCN) 2 to Cu[N(CN) 2 ]Br 
corresponds to increasing the electronic correlations or 
increasing of U/t. For example, it has been observed 
that as the pressure decreases m*/m e increases and the 
metal-insulator transition is approached 0-0]. Our cal- 
culated variation of the position of the dip with U is in 
qualitative agreement with the observed variation of the 
position of the dip with pressure (compare Fig. 4 in Ref. 
B). However, caution is in order because further exper- 
iments [ [l5[ find that the softening only occurs for waves 
propogating parallel to the layers when their polarisa- 
tion is perpendicular to the layers. Our model can only 
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explain this if modulation of the interlayer spacing has 
a much stronger coupling to the electronic charge den- 
sity within the dimers than modulation of the interdimer 
spacing. 

The family k-(BEDT-TTF) 2 X are particulary 
amenable to study the effect of electronic correlations 
on the molecular phonons because the dimer structure of 
the molecular crystal allows us to extract the electron- 
phonon coupling strength g from experimental data. The 
crystal structure is such that the BEDT-TTF molecules 
are arranged in pairs that are reasonably well separated 
from one another. There are three electrons per dimer. 
For each phonon mode on a molecule there is a sym- 
metric and an anti-symmetric combination on the dimer. 
By parity conservation, the anti-symmetric modes are 
infra-red active and the symmetric modes are Raman 
active. If there is a Holstcin model coupling for a single 
molecule the Hamiltonian for a dimer is lITrl 



Hp. 



t (4 



C2 



2 
i=l 



■ u)oa]ai 



where to is the transfer integral for hopping between 
molecules within the dimer. This is much larger than 
the hoppings t± and i 2 between dimers that appears in 
(0) H- The bonding and anti-bonding orbitals associ- 
ated with the dimer are separated by 2to, even in the 
presence of a Hubbard term. The Hamiltonian can be 
re- written as 



H ep — — to 



(c\c 2 



eld 



E 
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where a± = ^(o-i ±02) and n± — c\c% ± c 2 c 2 . The n + 
term is the same as the number operator, n^, appearing 
in (Q). It is the fluctuations in this term that produce 
the temperature dependence shown in Figures I and §. 
In contrast, the infrared active mode will not couple to 
these fluctuations and so should have no such tempera- 
ture dependence. 

If the number of electrons on the dimer is fixed then 
the symmetric mode has frequency lo . A second-order 
perturbative calculation JTJj predicts that for loq <C t , 
the infra-red frequency loi r will be smaller than the Ra- 
man frequency (u>r — wq), 



Zeyher and Zwicknagl [|T9| considered the frequency 
shift of phonons when they enter the superconducting 
phase. Within the framework of BCS theory, phonons 
with frequency lo q much smaller than the Fermi energy 
and much larger than the superconducting gap A will 
harden by an amount 



uj s - lu r _ 8g 2 N{0) 



L0 



UJQ 



(7) 



where N(0) is the density of states per spin at the Fermi 
energy. There should be no shift in the frequency of the 
infra-red active mode on entering the superconducting 
phase. Using A ~ 2fcsT c ~ 0.002 eV and the values 
for g in Table I we have evaluated the estimated shift 
in the phonon frequency for three different modes in 
Table | Eldridge et al §j| found that in k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br the u 60 (B 3g ) mode at 890 cm" 1 
hardened by about 0.2 % on entering the superconduct- 
ing phase. A shift in this mode was not detected in 
k-(BEDT-TTF) 2 Cu(SCN) 2 . No frequency shift in the 
other high-frequency modes was detected to within about 
0.1 %. However, one should be cautious about making 
quantitative comparisons between experiment and Table 
I since (0) is only valid for loq <C t and we have ujq ~ t. 
Pedron et al. did observe the hardening of modes with 
frequencies ranging from 27 to 134 cm -1 [ p0[ . For acous- 
tic phonons the softening will scale like (ujq/A) 2 and so 
be much smaller than the effects shown in Fig. |3[ 

In conclusion, we have shown how in a strongly corre- 
lated metal the redistribution of spectral weight over the 
scale of the band width with varying temperature can 
result in phonon frequencies with anomalous tempera- 
ture dependence. The effects involved are larger than 
those associated with the superconducting transition be- 
cause the latter only involves a redistribution of spectral 
weight over energies of the order of the superconducting 
gap which is much less than the band width. 
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Table I lists the details of the measured frequency shifts 
for three different modes for k- (BEDT-TTF) 2 Cu(SCN) 2 ■ 
For to — 0.3 eV [H these values are used together with 
(0) to evaluate the electron-phonon coupling g. The 
values obtained are in reasonable agreement with those 
found in a MNDO frozen phonon quantum chemistry cal- 
culation p| on a single BEDT-TTF molecule. 
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TABLE I. Parameters for three molecular phonons in k>(BEDT-TTF)2Cu(SCN)2. The experimental values are taken from 
Reference [^lj The difference between the infra-red frequency, loir and the Raman frequency, ljr, and equation Q is used to 
estimate (est.) the electron-phonon couplings g. These values are compared with values of g calculated by a quantum chemistry 
(QC) calculation Alu r (T) is the observed magnitude of the temperature dependence of the Raman-active mode between 
about 10 K and 300 K. Acusc is an estimate of the increase of the phonon frequency on entering the superconducting phase, 
based on equation (ph. 
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FIG. 1. Spectral densities of the metallic phase of the 
Hubbard model at half-filling calculated for the triangular lat- 
tice using dynamical-mean field theory. The bare density of 
states (U = 0) is shown as a dotted line. A peak at the Fermi 
energy due to a coherent band of Fermi liquid quasiparticles 
and two incoherent Hubbard bands situated at to — ±(7/2, are 
clearly resolved when the electron-electron interactions are 
larger than the metal bandwidth, W (U > W = 4.5*). The 
quasiparticle peak is strongly temperature dependent and its 
spectral weight becomes small as T > T* , T* , being the co- 
herence temperature. This strong temperature dependence of 
the spectral density is the origin of the temperature depen- 
dence of the phonon frequencies discussed in this paper. 
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FIG. 2. Temperature dependence of the real part of the po- 
larization for different frequencies and correlation strengths. 
Due to their coupling to electrons the phonons undergo a fre- 
quency shift that is proportional to this polarization. Signifi- 
cant temperature dependence occurs when two conditions are 
simultaneously satisfied: (i) u>o, is close to U/2 and (ii) elec- 
tron-electron interactions are sufficiently strong that there are 
well-defined Hubbard bands (compare Fig. 1). For U/t = 2 
the temperature dependence is negligible, as expected from 
the slight temperature dependence of the spectral density. 
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FIG. 3. Temperature dependence of the real part of the po- 
larization for u>o — 0.05* and different values of U /*. A dip ap- 
pears at the coherence temperature, T* , and it becomes more 
sharply defined as the electron-electron interactions increase. 
The position of the dip shifts towards lower temperatures as 
U/t increases. 
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